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,— , ■ Abstract 
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I . The nature of the gravitational interaction between ordinary and dark mat- 
ter is still open, and deviations from universality or Newtonian law may also 

I modify the standard assumption of collisionless dark matter. On the other 

I hand, obtaining a Yukawa-like large-distance modification of the gravitational 

' potential is a nontrivial problem, that has so far eluded a consistent realization 

■r^-j- ■ even at linearized level. We propose here a theory providing an Yukawa-like 

■ potential, by coupling non-derivatively the two metric fields related respectively 

. to the visible and dark matter sectors, in the context of massive gravity theories 

. where the local Lorentz invariance is broken by the different coexisting back- 

i grounds. This gives rise to the appropriate mass pattern in the gravitational 

Q I sector, producing a healthy theory with the Yukawa potential. Our results are 

^1 of a special relevance in the scenario of dark matter originated from the mirror 

' world, an exact duplicate of the ordinary particle sector. 
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1 Introduction 



The problem of obtaining a Yukawa-like potential in a consistent theory of gravity is a 
nontrivial task, unlike in (spontaneously broken) gauge theories, and attempts in this 
direction date back to 1939 when Fierz and Pauli (FP) added a mass term m to the free 
spin-2 action of the graviton [T]. However, the Lorentz-invariant massive FP theory 
is unfit to be a consistent modification of GR because of the van Dam-Veltmann- 
Zakharov (vDVZ) discontinuity [2J: even in the limit m — >■ the light bending is 
25% off than that is predicted from the GR and confirmed experimentally with an 
extremely high precision. Another theoretical problem is the fine-tuning needed to 
single out a ghost-free action at linearized level which however is probably spoiled by 
interactions and a sixth ghost-like mode starts to propagate making the whole theory 
unstable [3] and in any case unreliable below some (unacceptably large) distance scale. 
The problem was reexamined in the framework of effective field theory realising that 
the reason behind the misbehavior of FP massive gravity is strong coupling of the 
scalar sector [1]. 

It has been shown that the sickness of the FP theory has its roots in the Lorenz 
invariance [3]. Indeed, retaining only rotational invariance, one can avoid the vDVZ 
discontinuity and the propagation of ghost-like states P, [7j (for a different approach, 
see [8]). In these models gauge invariance is broken by Lorentz-breaking mass terms, 
so that the gauge modes that should start propagating, acquire a well behaved kinetic 
term, or do not propagate at all. What happens is that via Lorenz-breaking one can 
cure the 'spatial' problem of the discontinuity, while avoiding ghost-like propagating 
states. 

In the context of bigravity theories |9j, a suitable realization of the Lorenz-breaking 
(LB) massive phase of gravity can be obtained [TU]. In this approach, in addition to 
our metric field gi^^ coupled to the Lagrangian £i of normal matter (sector 1 in 
the following), one introduces another metric tensor g2fj,u related to a hidden sector 
2 (dark matter) with a Lagrangian €2- Therefore, the visible and dark components 
can be associated to separate gravities. The action of this theory consists of two 
Einstein-Hilbert terms and a mixed term V: 

S = Jd^x [v/^ {M^Ri + £1) + (M^ R2 + £2) + {gmf' V{X)j , (1) 

where Mi 2 are the "Planck" masses of the two sectors and e is some small mass 
scale which essentially will define the graviton mass trough see-saw type relation 
TTig ~ e'^/Mp, Mp being the Planck mass. The interaction potential V between the 
two metrics is assumed to be non-derivative and it can be always taken as a scalar 
function of Xj^ = gi°'g2au'^ the metric determinants are denoted gi and g2. The 
invariance under diffeomorphisms is not broken. Local Lorentz invariance, on the 
other hand, is spontaneously broken because in general there is no local Lorentz frame 
in which two metric tensors gi^u and g2f_iu are proportional. Nonetheless, because each 
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matter sector is minimally coupled to its own metric, the weak equivalence principle 
is respected and the breaking of local Lorentz invariance is transmitted only through 
the gravitational interactions. Once a flat rotationally invariant (double) background 
is found, a Lorentz-non invariant mass term for the gravitational perturbations arises 
in a natural way, by expanding the total action ([T]) in the weak field limit [TOj . 

We point out that when the interacting potential V is absent the (gauge) symme- 
try is enlarged: one can transform (71^,^ and g2^u by using two independent diffeomor- 
phisms. The interaction potential V leaves unbroken the common diffeomorphisms 
group, corresponding to general covariance, and hence the gravitational fields always 
include a massless sector including normal graviton. In the massive sector on the 
other hand one finds a massive graviton, and when Lorentz-invariance is broken there 
are no additional propagating modes: in particular vector and scalar degrees of free- 
dom do not propagate [TU] . The Newtonian potential is modified in the infrared, but 
it is not Yukawa-like. In fact, at linearized level, the deviation from a 1/r potential 
is a linearly growing term [101 Cj ^ 

In order to have a massive phase with a Yukawa-like potential, bigravity must then 
be enlarged. In this paper we generalize the above construction and show that one can 
use a further rank-2 field g^^i, as a Higgs field to achieve the Yukawa potential. The size 
of the fluctuations of gsf^u is controlled by the 3rd " Planck" mass M3 entering in its EH 
action. We will show that, in the limit M3 ^ Mi^2, fl's/^i/ can be consistently decoupled 
and one is left with an effective bigravity theory with a Yukawa-like component of 
the gravitational potential. The tensor (^3^,^ plays the role of a symmetry-breaking 
field, communicating the breaking of Lorentz invariance to gi^j^i, and g2^u and thus 
introducing Lorentz-breaking mass terms to their fluctuations. Even in the limit of 
the decoupling of gs the resulting phase of gravity features a Yukawa-modified static 
potential while still avoiding any propagation of ghosts and the vDVZ discontinuity. 

This situation can open new possibilities for the nature of dark matter. In the 
present paradigm the visible matter amounts only to about 4% of the present energy 
density of the Universe while the fraction of dark matter is about 5 times bigger. 
Cosmological observations are consistent with the hypothesis of cold dark matter. 
On the other hand, the flattening of galactic rotational curves can be also explained 
by the presence of cold (coUisionless) dark matter distributed, differently from the 
visible matter, along the galactic halos. The implicit assumption behind this scenario 
is that gravitational interaction between the two kinds of matter is universal, and 
that it is Newtonian. Relaxing these hypotheses may radically modify our view and 
phenomenological modelling. 

One of the intriguing possibilities is to consider dark matter as a matter of a hidden 
gauge sector which is an exact copy of the ordinary particle sector, so that along with 

""^This term breaks perturbativity at some large distance r > riR, but remarkably this behavior 
is cured by the non-perturbative treatment [TT] . 
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the ordinary matter: electrons, nucleons, etc. the Universe contains also the mirror 
matter as mirror electrons, mirror nucleons, etc. with exactly the same mass spectrum 
and interaction properties. Such a parallel sector, dubbed as mirror world ^12j, can 
have many interesting phenomenological and cosmological implications (for reviews, 
see [IH]). In particular, the baryon asymmetry in both sectors can be generated via 
the out-of-equilibrium, B — L and CP violating processes between the ordinary and 
mirror particles [H] which mechanism could naturally explain the proportion between 
the visible and dark matter fractions in the Universe. Such processes can be induced 
by the some very weak interactions between the ordinary and mirror fields that on the 
other hand can induce the mixing terms between the neutral particles of two sectors, 
as e.g. kinetic mixing for photons [T6j or mass mixing in the case of the neutrinos 
and neutrons [T7] . 

Mirror matter, dark in terms of ordinary photons and coupled with ordinary 
matter via common gravity, can be a viable candidate for dark matter. As it was 
shown in [15j, the cosmological observations on the large scale structure and CMB 
are consistent with the mirror dark matter picture. However, the essential problem 
emerges at the galaxy scales. It is difficult to understand how the mirror matter, being 
as collisional and dissipative as normal matter, could produce extended galactic halos 
and thus explain the galactic rotational curves. 

In this paper we show that the new possibilities can emerge if the mirror sym- 
metry is extended also to the gravitational sector in the form of the action ([T]), the 
normal and mirror matters having separate gravities related respectively to the met- 
ric fields Qi^ij and g2^v while the Lorentz breaking is induced by the third dynamical 
metric (73^,^ with its "Planck" mass M3 much larger than the ordinary Planck mass 
Mp. Two gravities - one massless and another with a nonzero mass m, lead to the 
Yukawa-modified gravitational potential along with the normal Newtonian term. The 
potential felt by the probe particle of the type 1 (normal matter) at the distance r 
from the source is 

G 

(p{r) = — [(mi + ma) + (mi - ms) e"'^/'"'"] , (2) 

where G is the Newton constant, mi and m2 are respectively the masses of the visible 
(type 1) and mirror (type 2) matter sources, and = m~^ is a Yukawa length scale. 
Hence, at small distances, r <^ r^, the gravitational forces are not universal between 
two sectors: the normal and mirror matter do not see each other. At distance r <^ 
a normal test mass interacts only with mi through the ordinary Newton potential. 
But at large distances r ^ Vm the gravity becomes universal, and test particle feels 
both ordinary and dark matter sources (mi -|- ma) with an effectively halved Newton 
constant G/2. The main result of this work is to reproduce the potential ([2]) in a 
consistent model of gravity. 

This scenario can have interesting astrophysical implications. One can show [23] 
that it allows to reproduce the galactic rotational curves even if dark mirror matter 
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has the similar "clumped" distribution as the normal matter, as it is expected from 
its dissipative character. 

The paper is organized as follows: in section 2 we review the linearized analysis 
of bigravity theories, to be used as building blocks for the model, and describe their 
Lorentz-breaking and Lorentz-invariant phases. In section 3 we describe the model 
and show how a Yukawa potential arises in the limit when the additional metric is 
decoupled. In section 4 we discuss the findings. Finally, appendices A and B contain 
the detailed expressions for the graviton mass matrices, the details for a specific 
interaction potential and the general expression of the Yukawa-like potential. 



2 Bigravity: A review of the Linearized Analysis 

In for bigravity generically one can find bi-flat SO {3) preserving vacuum solutions \TU\ : 

^1/.,. = ^/.i. = diag(-l, 1,1,1) 
g2^„y = fif,u = t^Miag(-c^ 1, 1, 1); 

we have set the speed of light in our world (sector 1) to be one in natural units, 
whereas c is the speed of light in the hidden sector 2 and a; is a relative constant 
conformal factor. Once V is given, c and u can be computed by solving the equations 
of motion following from ([T]), and if c 7^ 1, Lorentz symmetry is broken. Consider the 
linearized theory obtained by expanding the total action ([T]) at quadratic level in the 
metric perturbations around the bi-flat background 0: 

gifiu = Vt^u + hi^y , g2f,u = Vfiu + t^^^2Mv • (4) 

The gravitational perturbations /ii^jy and h2fj,u interact with matter 1 and 2 through 
their conserved EMTs, respectively T^^ and T2'' . Since the background preserves 
rotations, it is convenient to decompose the perturbations hafj.u [a = 1, 2) according 
to irreducible 5*0(3) representations 

^aOO = "^a 5 h,aOi = Uai + diVa , , , 

haij = Xaij + diSaj + OjSai + Oidjaa + 6ij Ta . 

For each perturbation one has a gauge invariant transverse traceless tensor Xaij, two 
vectors and four scalars. The quadratic Lagrangian C reads 



•£ -^kin ~l~ -^mass ~l~ -^src 5 (6) 

+ (j)'JCA T - -t'IC {C^A - 3d^) T. (8) 
2 



C^in = -x%nC^^-dl)x:j-7^w\lCAw,+ (7) 
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We have introduced a vector notation for the fluctuations: h^^ = {hi^^, /i2/ii/)*, Xij = 
(Xiij,X2ij)* and the following 2x2 matrices: C = diag(l,c), /C = diag(l,K) and 
K = M2/MIUJC. Also, in the kinetic term, coming from the expansion of EH terms, 
the fluctuations enter only through the gauge invariant combinations Wi = Ui — dtSi, 
(j) = ^ - 2dtv + d^a. The mass term £mass is produced by the expansion of the 
interaction potential V. Finally, Csrc describes the gravitational coupling to conserved 
sources associated to matter fields: 

A,, = C-'W, - rj-^cl> - Ti'^T - Tl^x^, . (9) 

Clearly £kin and £src are gauge invariant. For the bi flat background ([3]) the mass 
term Lm has the following form 

>Cmass = J [hlo-^ohm + 2hliMihoi - h\jM2hj + h^M^hu - 2/1*^X4^00] (10) 

and the explicit value of the mass matrices can be easily computed for any given V . 

It is however crucial to realize that due to linearized gauge invariance the mass 
matrices have the following property jlO] 

-^0,1,4^2) =0' -Mi,2,3Q=0, A^*Q=0. (11) 

Thus general covariance forces the mass matrices to be at most of rank one. 

Lorenz-Invariant (LI) phase. In this case an FP graviton mediates Yukawa-like 
potential. Indeed, when c = 1, two conditions in ffTTl) coincides , allowing a non-zero 
M.I and all mass matrices are rank one and proportional: 

Mo = \oV, Mi = M2 = \2V, M3 = M4 = i\2 + \o)V , V = 

(12) 

After introducing a canonically normalized graviton field h^^^ = JC^^'^h, the mass 
matrices can be diagonalized by a rotation of an ang le with tan 7? = k^/^ ^ M2/0JM1, 
leading to a massless and massive graviton eigenstates. This latter has a standard 
Lorentz-Invariant mass term and to avoid ghosts one has to choose Aq = 0, leading 
to a Pauli-Fierz mass term. Then, the massless graviton interacts with both matter 
sectors and effectively acts like standard Newtonian gravity in the weak field limit. 
The massive graviton on the other hand is Yukawa-like and thus modifies the static 
potential at scales larger than m~^, where m = e^A2^^| sin??|/Mi is the graviton mass. 

In the most interesting (mirror) case, when Mi = M2 = M, u = 1 and so 
tan-i? = 1, we obtain that a static potential for a test particle of type 1, generated 
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by the point-like sources of mass mi (type-1) and m2 (type-2) at the same point, is 
given by 

01matter(^) 



4 

(mi + 1712) + - e (mi - 1712) 
o 



(13) 



327rM2 

Therefore, the presence of the massive graviton state mediating Yukawa-hke terms 
makes the effective Newton constant distance dependent: the Newton constant mea- 
sured experimentally via the gravitational interaction between the type-1 test bodies 
at small distances r ^ m~^ should be identified as G = Guv = T/dGirM"^ while 
at large distances r ^ m~^ it effectively becomes G//? = l/327rM2 = 3G/7 and is 
universal between the type-1 and type-2 matters. On the other hand, at r -C m^^ 
the gravitational forces between the type-1 and type-2 bodies are repulsive, with 
Gjjy = —G/7, which in itself is indication of the instability of the theory. However, 
more serious problem is related to the vDVZ discontinuity. The static potential felt 
by the photon is 

<^iiight(r) = [(mi + ma) + e"'"^ (mi - ma)] . (14) 

Therefore, for the light bending at distances r <^ we have Gught = l/lGvrM^, 
and thus G/Gught = 7/6. This discrepancy is somewhat milder than in the FP 
theory where we have G/Gught = 4/3; anyway the deviation from the GR prediction 
G/Giight = 1 is unacceptably large and it is clearly excluded by the post-Newtonian 
gravity tests [18] . 

The problems can be softened if the two sectors are not symmetric. Mi 7^ M2 
and the mixing angle between two gravities is enough small. In this case the static 
potentials respectively for the test body and test photon of the type 1 read: 

(pidr) = ^Q^j^2 [{mi + m2)+^e-"''' {m^ tan^^^-ms)] , (15) 

where ^ = 4/3 for a test body and ^ = 1 for light. Therefore, at small distance 
(r <^ m~^) the Newton "constant" is Guv = G{1 + i/Staii^ -d), at large distance 
(r ^ m~^) it tends to G. 

The ratio of (JT5l) at small distances defines the post-newtonian parameter 6: 



lim 



= 4/3) 

0i?(e = i) 



= l + ^sm'^^. (16) 

m2=0 3 



The current light bending experiments put a constraint 6 = 1.0000 ± 0.0001 [18], and 
for GR (5 = 1, so that the limit of vanishing graviton mass reveals the well known 
vDVZ discontinuity |2] of Pauli-Fierz massive gravity. In our case the mixing angle 
i) controls the size of the discontinuity. 

When M2 ^ Ml, we have 7r/2 and 6 = 4/3, unacceptably large. In this limit 
sector 2 is very weakly coupled, and the discontinuity is mainly shifted to sector 1, 
that approaches a normal Fierz-Pauli massive gravity. 
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Conversely when M2 -C Mi we have ^9 — and the discontinuity is shifted to 
sector 2; /i+ and almost coincide with hi and /12 and gravity is stronger in sector 2. 
The experimental bound on S translates into ?9 ~ 0.02, that amounts to roughly M2 — 
dMi. In this case, if m2 is interpreted as dark matter, it gives a sizable contribution, 
increasing the gravitational force in the region r > mT^ . Notice incidentally that 
for small r, the potential is repulsive. This result contradicts observations in the 
gravitationally bounded systems as cluster and galaxies, for this reason is ruled out. 

Lorenz-Breaking (LB) phase. In this phase, 07^!, conditions ( |TT|) imply that 
M.i=Q and for other masses one has 

Mo = \^C~^VC-\ A^2,3 = A2,3^, M^ = KVC-\ (17) 

In this situation all the scalar and vector perturbations become non-dynamical [TU] . 
The vanishing oi M.i in the LB phase is the reason why no ghosts or tachyons appear 
in the theory and only gravitational waves propagate. However, this is also the reason 
behind the absence of Yukawa-like gravitational potential. The resulting modification 
was studied in detail in [Tlj both at linear and non-linear level. 

3 Three metrics: Effective Higgs Phase 

In order to find a phenomenologically healthy Yukawa phase, we introduce one more 
rank-2 field which couples with both gi and g-^ 

S = jd^x [M^Ri + £1) + (M^ i?2 + C2) + Ml V^i?3 

(fi'i5'2fi'3)^^^^(5'i,fi'2,5'3) • (18) 

The only non trivial tensors that can be formed are: X12 = gi^g2, ^is = 91^93, 
X23 = 92^93: that satisfy the identity X13 = X12X23. Therefore V can be taken as a 
scalar function of two of them. 

We also introduce in ( |T8|) a discrete symmetry under the exchange 1^2, so that 
the potential V is symmetric and the two sectors 1, 2 have equal Planck masses 
M. The third Planck mass on the other hand will be eventually taken much larger, 
M3 ^ M, and the fluctuations of the third field will be effectively decoupled. 

The first step is to find a suitable background. As for bigravity, we look for fiat 
solutions, for which a consistent LB ansatz is the following 

9if,u = g2t,u = Vf^u = diag(-l, 1, 1,1) 

(19) 

g3^lu = Vfiu = ^'^ diag (-c^ 1, 1, l) , 



^For the simplicity, we assume that the third auxihary sector is purely gravitational and does 
not contain the respective matter. In principle, any tensor condensate e.g. emerging via a strongly 
coupled hidden gauge sector can be also used for inducing the Lorentz-breaking background [H] . 
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so that Xi2 = I and X13 = rj ^7). The background ( |T9l) is a solution of the equation 
of motion if the following equations are satisfied: 

dV dV 
yi + 6-— X21 + 6-— X31 = 
5X21 dX^i 

dV dV 
V"I + 6- X12 + 6- X32 = 

dV dV 
yi + e-— X13 + 6-— X23 = 0, (20) 

CA13 0'A23 

— 1 



where X^^ = X^^^ . Then, due to the identity 



dV dV 

-Xab = —777^ ^fea , (21) 



and due to the 1^2 exchange symmetry of the EOM and of the background, we have 
dV/dXi2 = and the EOM reduce to 

dV 

V^O, — = 0, (22) 

These are three independent equations for the two parameters u and c: thus one 
fine-tuning is needed to have the present fiat solution. This fine tuning is analogous 
to the cosmological constant in standard GR, and can be easily realized for instance 
by introducing a cosmological constant in sector 3. 

Once a background solution is found, one can study small fluctuations defined by 

gii^iu = f]nv + hi^y , g2^v = rj^y + h2^v , Q^^iv = Vfiu + t^^^s^i^ • (23) 

The structure of the quadratic Lagrangian for the fluctuations is the same as in ([6])- 
([9]) except that now the tensor, vector, scalar and source fields all have 3 components, 
h2i,v,h'if,^). Also, 

K. = diag(M2, M^, M^/u'^c) , C = diag(l, 1, c) (24) 

and the masses Aii are 3x3 matrices, entering the usual mass Lagrangian: 

jCmass = hl^Mohoo + 2hliMihoi - h\.M2hij + hl^Msha - 2h\^MihoQ . (25) 

Diagonal diffeomorphisms invariance constrains the form of these matrices: 

A^o,i,4 I 1 1=0. (26) 
c2 
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From these conditions and from the 1^2 symmetry it follows that the matrices can 
be written as the following combinations of projectors 



where 

Vl2 





Mo 


= aoVi2 + boC-\Vi3 


+ 7^23)^-' 








Ml 


= aiVi2 










M2 


= a2Pl2 + b2{Vl3+V23) 








M3 


= 03^12 + b3{Vl3 + V23) 








Ma 


= 04^12 + 64 (7^13 + V23) 






( 1 


-1 q\ 






( 1 -l\ 




{ 


\ 


-1 


1 




Vl3 = 





P23H 


1 


-1 


I 


y 






I -1 1 J 




VO -1 


1 ) 



(27) 



(28) 



and Oj, bi are constants that depend on the quadratic expansion of the interaction 
term V. 

Since we are interested in the gravitational potential we will first focus on the 
scalar sector. The full Lagrangian is: 



scalars 



At ^2 
1 

+ - 



■^Mo'^ - lAvMxv - (r + A(7)A^2(t + Act) - 2tM2T 



+(3r + Aa)A<3(3T + Aa) - 2(3t + l\a)M4} 



00 



(29) 



In order to disentangle the different fluctuations we decompose the system by defining 
a 'tilded' basis where the 1,2 fluctuations are rotated: 



V 



[i^,v,a,T]^S^,v,a,f^, Mi^S'MiS S 
In this basis the mass matrices take the block-diagonal form 

Ml 

Ma^ 



^ ^ 

4s 



(30) 



V2 / 





' 2ao + 60 





Mo^ 1 


bo ■ 


-b,/c' 




. -bo/c^ 


bo/d" 




/ 2^2,3 + 62,3 







62,3 


-62,3 




V -62,3 


62,3 




'204 + 64 

64 -bi/c^ I (31) 
-64 64/c^ 
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Because the kinetic structure commutes with S, we see that in the new basis 
the system sphts into two sectors: a single massive gravity plus a bigravity sector, 
associated with the 2x2 sub-matrices in (IHTI) . Due to the third background g^, in 
both sectors the mass pattern is Lorentz-breaking. The first sector can be analysed 
as in [6] , while for the second the analysis of applies. This allows for a consistent 
theory, free of ghosts and of instabilities at linearized level. 

Indeed, in the single massive gravity sector ghosts can be avoided, if the relevant 
entry 1-1 in Mq vanishes. We have thus the condition: 

ao = -bo/2 . (32) 

The bigravity sector on the other hand is automatically free of ghosts as shown in [TU] 
thanks to the vanishing of A^i in the relevant block. 

At this point we can study in the new basis the static gravitational potential 
associated in each sector with the gauge invariant field 4>a = i^a — 2dtVa + d^da 
(a = 1,2,3). It is convenient to define also the rotated and M^-normalized sources 
t^. = St,, = S{T,JM^). 

The field 0i is separated from the bigravity sector and gives the Yukawa-like static 
potential. It turns out that in general 0i is a combination of two Yukawa potentials, 
with two parametrically different mass scales (see appendix |A] for the details). For 
simplicity, by tuning the parameters one can also have a single mass scale: 

ti 

with = 3(2a4 + 64)-— ^ . (33) 



In this sector, in addition to the propagating massive graviton (two polarizations) also 
a vector and a scalar field propagate (respectively two and one degrees of freedom). 
All these fields are massive with mass given by the relative 1-1 entry of M.2- The 
vector and the scalar can have well behaved properties, i.e. no ghosts when condition 
is enforced. In [0] it was also argued that the scale of strong coupling is high 



enough, coinciding with A2 — V Mm, with m the characteristic mass scale in this 
sector. 

For the remaining bigravity sector the gravitational potential can be computed by 
solving the equations of motion as in [10]. The result is 



i^200 + ^2m , 2 '''200 



where 

.2 



■ 3bl + bo{b2-3h) 
,'2ibl + boib2-h))\ 



(36) 
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When M3 3> M, the third sector has a sub-leading impact on the other gravitational 
potentials. In the limit M3 — 00, the third sector decouples and just produces a 
LB fixed background fj. Going back to the original basis, the potentials are: 



^100 + tin + ^^200 + t 



''2m 



4oo 



4A 

tin + ^200 



+ 



tlOO + t 



Hi 



200 



t 



2ii 



''2m 



'100 



4A- 

~ tin 



- ^200 



+ fi 



2^100 + t 



200 



2A2 



t2ii 



-/i 



M 
Ms 



4A 

2 



4A - 2m2 



2^100 



^^200 



2A2 



(37) 



CUJ^{t 



100 



''200 J 



A2 



The final potentials 0i_2 contain a Newtonian term, a Yukawa-like term, and a linearly 
growing term, originating from /i^/A^. 

This latter linear term is the same appearing in the bigravity found 
in [ini 120] • It would invalidate perturbation theory at distances larger than rj^ ~ 
Gn^Mj, from a source M^, [TU], but remarkably the full nonlinear solutions found 
in [H] shows that its linear growth is replaced by a non- analytic power ~ r'^, where 7 
depends on the coupling constants in the potential. Moreover, in the full solution for 
a realistic star, also the magnitude of this new term is proportional to /i^, therefore 
the effect can be eliminated by setting = 0. This can be achieved by simple fine- 
tuning, or by adopting a particular scaling symmetry of the potential, as discussed 
in [To]. We can thus obtain a pure Yukawa modification of the gravitational potential, 
by setting /i^ = 0, that here amounts to the condition bo = —364/(62 — 363). 

The analysis of vector modes is identical to that carried out in [6J for the single 
gravity sector and to the one of [TU] for the bigravity one. In the single-gravity sector 
there is a vector state propagating with nonlinear dispersion relation: at high energy 
its speed is (2a2 + 62)/(2ai -|- bi) and at low momentum it has a mass gap given by 
62/^^. In the bigravity sector on the other hand vector states do not propagate. 

The analysis of tensor modes is similar and is best carried out in the original basis. 
In the limit of M3 —>■ 00, the equation of motion for the canonically normalized fields 
becomes: 















1 



















'62 + 02 62 - 0,2 
62 - 02 62 + O2 





Xlij\ 




(tlij 


X2ij 




t2ij 


Xsij/ 







{31 



where □ = rj^'^d^d^, □ = fj^^d^dy and we used the form of the projectors (j271l . We 
see that the massless spin two state decouples (it is a superimposition of mostly X3) 
and we are left with two massive gravitons, with two polarizations each, travelling at 
the normal speed of light. Their mass matrix can be diagonalized, and the resulting 
graviton masses are m?g^ = (2a2 + 62)^'^/^^, = b2e^ /M"^. 
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As an explicit example, consider the simplest case of a potential quadratic in 
Xi2, Xi3, X23 plus two cosmological terms, satisfying the 1^2 exchange symmetry 
(taking, for simplicity u = 1): 

Vigi,g2, gs) = + 6 (Tr [X^g] + Tr[X|3]) + 6 Tr[Xi3X23] + 
^3 (Tr[Xi3]' + Tr[X23]') + ^4 Tr[Xi3]Tr[X23] + 

{Tt[Xu]' + TriXf^T) + ^6 {Tr[Xl,] + Tr[(Xr/)2]) + 
er ((detXi2)-^/^(detXi3)-^/^ + (detXi2)^/^(detX23)-^/^) + 
e8(detXi3)^/^(detX23)'/' (39) 

The EOM for fiat backgrounds require to solve for three constants (e.g. .^3, 1^7, ^g) and 
then the coefficients of the projectors in the mass matrices aj's and 6j are a function 
of the remaining coupling constants (see appendix [B|) . 

The no-ghost condition bo = — 2ao, the condition for the absence of the linear term 
/i^ = and the condition for having a single Yukawa scale (see Appendix |A]) , can be 
solved for ^1,2,4 and reduce finally the dependence to only four couplings. Then, the 
Yukawa scale m and the graviton masses m'^^ and m^2) oiily depend on ^5 and ^q: 

m'^ = [Po(c)^5 + go(5)^6]^ 



^ii = [Pi{c)^5 + qi{c)^6]j^ (40) 



^g2 = [P2(c)6 + g2(c)^6] 



where pj(c)'s and gj(c)'s are given in appendix [BlFI 

To summarize, in the limit where the third metric is decoupled the theory has two 
massive gravitons and the potential felt by a test particle of type 1 is: 

'^iW = — ^ + — ' (41) 



where G = l/lGvrM^. This shows that the vDVZ discontinuity is absent, and we 
have obtained the potential ([2]) while avoiding the troubles of the Lorentz-invariant 
FP theory. 



4 Conclusions 

Motivated by the interesting possibility to relax the assumption that dark and visible 
matter feel the same gravitational interaction, in this work we have addressed the 

■^When ^6 = and all the masses above depend only on ^5, one can check that they are positive, 
for 1.41 <c< 2.05. 
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possibility to obtain a Yukawa-like large-distance modification of the gravitational 
potential, while avoiding ghosts or classical instabilities. 

The request to generate an Yukawa potential from a consistent theory led us to 
consider Lorenz-Breaking backgrounds in enlarged models of bigravity. For example, 
bigravity theories while giving rise to a healthy Lorentz-Breaking massive phase, do 
not produce a Yukawa potential. Here we have generalized this picture and have 
shown that if an additional field gs is introduced, a Yukawa modification is allowed. 
The extra field can be harmlessly decoupled by freezing it to a Lorentz violating 
background configuration. The two remaining sectors represent two interacting mas- 
sive gravities, of which one features a Yukawa potential. This pattern then leads to 
the desired modified gravity, because standard matter (type 1) couples to all the mass 
eigenstates. 

On the technical side, the price to be paid to solve the problem is that two fine- 
tunings are needed, one to have a ghost-free spectrum at linear level, the other to 
avoid the linearly growing potential. The first one has be shown to follow (in single 
massive gravity theories) from extra unbroken partial diffeomorphisms invariance [21] 
and it would be interesting to extend that symmetry arguments also to the present 
model. The other can also be understood as the consequence of a scaling symmetry 
of the potential [TO] . 

Let us note also that while the theory presented of three rank-two fields only 
propagates 9 well behaved modes at quadratic level (three spin-2 with 2 polarizations 
each, one spin-1 with 2 polarizations and one scalar) one may expect that the total 
number of propagating modes would be 18, i.e. 3x10 minus 8 by unbroken gauge 
conditions minus 4+4 for the broken relative diffeomorphisms (a la Proca). The 
missing 9 modes could then propagate at non-linear level. The real non-perturbative 
question, to be addressed in future work, is then at which scale these non-linearities 
would show up. 

The resulting setup, featuring two separate metric fields responsible for gravity 
for the visible and dark matter, allows to consider also coUisional and dissipative 
dark matter, as mirror matter, if the potential felt by ordinary matter, and generated 
by the ordinary and dark matter sources of mass mi, m2, is distance dependent as 
m The result is very different from a standard picture when normal and dark 

matters both have an universal Newtonian gravity; in fact the potential fl4ip can be 
used to fit galactic rotational curves using similar density profiles for the visible and 
dark sectors, alleviating the problems of profile formation [23] . It interesting to note 

* Let us remark also that the weak equivalence principle does not exclude the possibility of direct 
(non-gravitational) interactions between the normal (type 1) and dark (type 2) matter components. 
To the action ([T]), besides the mixed gravitational term V, the mixed matter term J (i'*x((/i(72)^^'*'Cmix 
can be added with the Lagrangian >Cmix including for example, the photon kinetic mixing term 
eF^'^ ¥2^1, [IH] or the neutrino interaction terms [T7]. This also makes possible the direct detection 
of dark matter via such interactions [22] . with interesting implications. 
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that the effective Newton constant relative to the type 1 - type 1 and type 2 - type 2 
interactions is distance dependent: GAr(r -C m~^) = G and GAr(r 3> m~^) = G/2. 
Let us remark also, that as far as at large cosmological distances only the massless 
gravity is effective with a halved Newton constant G/2, the observed Hubble constant 
would imply for the total energy density of the universe twice as bigger than in the 
standard cosmology when the Newton constant at the cosmological distances remains 
the canonical G, i.e. now we must have pcr = 3H^/AnG instead of pcr = 3i7o/87rG 
impled by the standard cosmology. 

We then conclude that at linear level a Yukawa modification of the Newtonian 
gravitational potential is possible and it also opens up to the possibility to have coUi- 
sional dark matter, coupled to ordinary matter via a modified gravitational interaction 
in a physically nontrivial way. 
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A General Yukawa-like potential 

The degrees of freedom in the single gravity sector consists of a metric fluctuation 
with mass term that we can write as 

^mass = {^IhmhoQ + 2mlhoihoi - mlhijhij + mahuhjj - 2m\hoQhii) . (42) 

In our case effectively mo = 0. In this case, if nii ^ 0, there is a (healthy) propagating 
scalar degree of freedom (r) as well as a healthy propagating vector [6j. The scalar 
perturbations obey the equations: 

2Ar - m2(Aa + 3r) = too (43) 
2doT-mlv = ^dotoo (44) 

2dQr — mlAa — m\Tm\Aa + ?,m\T — ra\{(t) + 2dQV — d^a) = -^c^otoo (45) 

3 

2A0 - 2Ar + 2mlAa = tu - — 9^too ,(46) 

where t^i, = T^y/M"^. These can be solved with respect to to get the static Newto- 
nian potential. One finds 

^ _ (too + (Cl - 1)C2A + [tu + too(3Cl - l)C2]C2ml ^^^^ 



2(Ci - 1)C2A2 + (4C2 - l)m2A - 3C2ml 
15 



5 



with (i = m\/m\ and C2 = Tn^/m?^. The potential can be spht in two Yukawa-hke 
terms: 



2A - ml 2A - m 



where 



t± = i(too + tn) (^1±^^ ±too (^5-0 , (49) 

^1 = rnl , 5= Vl + 8C2(3CiC2-C2-l). (50) 

Recall [6] that the conditions (2 > 1/4 and Ci < 1 ensure that the theory as no 
derivative instabilities neither in the UV nor in the JR. Moreover, if (i > (8C| + 8(2 — 
1)/24C|, the two masses m± are real and positive, and the theory has no instabilities 
also at intermediate scales. Accordingly, the potential is the sum of two "genuine" 
Yukawa-like terms. 

Finally, if Ci = (1 + C2)/3C2, then 6 = 1 and t_ vanishes, so that one is left with a 
single Yukawa potential: 

too + tii 



2A - 3ml 



(51) 



B Explicit solution for potential (1391) 



Mass coefficients as a function of the coupling constants for the potential (l39l) . after 
solving the EOM 



ao 
bo 



.^-^-^-^(9c^ + 2)6 + ^ + ^ 
2 9 72 18 ^ "^ 9 9 

eo cM39 - 23c^)ei c^(39 - 23c^)6 ^5 ^6 

"72 18(c2 + 3) 36(c2 + 3) 9 9 



ai = 4^5 

a2 = ^-^5 

b2 = (c^-i)ei + ^(c^-i)6 

= —9 ^-72-2+^ + ^ ^^2) 

6 ^ (5c2-6c)ei ^ (5c2 - 6) 6 6 46 
03= - — H tt; H 



04 

&4 



72 18 36 9 9 

26c2 6c' 6c' 6 6 , 506 



9 9 2 72 9 9 

6 (13c^ + 3c2)6 (13c^ + 3c2)6 6 46 
'72 18(c2 + 3) 36(c2 + 3) 9 9 
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Functions p appearing in the graviton masses (HOl) : 

= -1 [6 (c^ ( (3 (I850c^ - 7725c^ - 31099c^ + 1545070^ - 92547) + 7Ci- 

168318) - 5Ci + 117936) - I8C1)] (53) 

go = 7^ [12 (c^ ((3610c^° - 15312c^ - 58045c^ + 296415c^ - 1874610^+ 

UCi - 115371) - lOCi + 88452) - 36Ci)] (54) 

Pi = --^ [2 {95c^^ - 56740^^ + 21090c^^ + 950530^^ - 447746c^°+ 

243567c^ + 194157c^ - 37Ci (6c* - 125117c^ + 5c^ + 42))] (55) 

qi = [2 (380c^* - 7716c^^ + 22284c^^ + 1231140^^ - 529730c^° + 334098c^ - 148+ 

Ci (24c* - 205631c^ + 200^ + 168238626))] (56) 

(95ci° - 324c* - 405c6 + 378c^ + QCi) C4 , , 

P2 = (57) 

4 (95cio - 324c* - 405c6 + 378c^ + 6^1) C4 
12 = ^ , (58) 



where 



Ci = c^ (5c^ - 26c2 + 21) Vl3c4 + 42c2 + 9 (59) 

C2 = (1404- 1773c2 + 6c^ + 107c^)(c^(-2 + c^)(-21 + 5c2)) (60) 

C3 = c*(58968 - 117990c2 + 62235c^ - 4557c^ - 3287c* + 535c^°) (61) 

C4 = 2(c^ + 2c2-3) (62) 

C5 = c''(-29484 + 44253c2-8991c^-2217c^ + 535c*). (63) 
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